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Study of Magnetic Excitation in Singlet-Ground-State Magnets 
CsFeCls and RbFeCla by Nuclear Magnetic Relaixation 
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The temperature dependences of spin-lattice relaxation time Ti of ^^^Cs in CsFeCls and *^Rb 
in RbFeCls were measured in the temperature range between 1.5 K and 22 K, at various fields 
up to 7 T applied parallel (or perpendicular) to the c-axis, and the analysis was made on the 
basis of the DCEFA. The mechanism of the nuclear magnetic relaxation is interpreted in terms 
of the magnetic fluctuations which are characterized by the singlet ground state system. In the 
field region where the phase transition occurs, T^^ exhibited the tendency of divergence near 
Tn, and this feature was ascribed to the transverse spin fluctuation associated with the mode 
softening at the /('-point. It was found that the damping constant of the soft mode is remarkably 
affected by the occurrence of the magnetic ordering at lower temperature, and increases largely 
in the field region where the phase transition occurs. 
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§1. Introduction 

Hexagonal compounds CsFeCls and RbFeCla are typ- 
ical singlet ground state magnets. The magnetic system 
of Fe^+ ion which is described in terms of a fictitious spin 
of S" = 1, is characterized by the singlet ground state and 
the excited doublet states due to a large positive single- 
ion anisotropy along the crystallographic c-axis. Here- 
after, these two compounds are abbreviated as CFC and 
RFC, respectively. Because of the competition between 
this anisotropy and the dominant ferromagnetic intra- 
chain exchange interaction along the c-axis, these com- 
pounds exhibit various magnetic properties. The com- 
mon features are the softening of magnetic excitations, 
which gives rise to 3-dimensional long-range order (3D- 
LRO) such as the commensurate and incommensurate 
phases. On the other hand, there appears different fea- 
tures originating from a difference in the relative mag- 
nitude of the anisotropy with respect to the intra-chain 
exchange interaction, even if it is small. 

So far, magnetic properties such as an excitation spec- 
trum including the characteristic softening effect and a 
phase transition, have been studied extensively by the 
measurements of magnetization, ^"^^ specific heat^'^-* and 
inelastic neutron scattering. The phase diagrams 
have been determined by the specific heat measurements 
for both H//c^^ and _ff_Lc.^' In the case of CFC, the sys- 
tem remains non-magnetic down to zero K at zero field. 
While, when the external field is applied along the c-axis, 
one of the excited doublet decreases with increasing field 
and the energy levels of the ground state and the excited 
state cross each other at the field He (=7.5 T). Then the 
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field induced ordered state appears in the field range of 
4 T<i?<ll T around In the case of RFC, the 

magnetic phase transition occurs at zero field at Tn ~2.6 
K, and the ordered state disappears when the value of 
the external field is beyond 2.2 T in case of HJ-c}'^''^'^'^'' 

In particular, it should be noted that the softening of 
magnetic excitation was observed in the inelastic neu- 
tron scattering experiments in RFC.^-* It was recognized 
that the phase transition is driven by a specific mode 
of magnetic excitation which decreases to zero at the 
phase transition point (Tn) with the ordering vector q = 
(1/3, 1/3, 0) (if -point). Suzuki developed the dynamical 
correlated-effective-field approximation (DCEFA) fo- 
cusing on the problem of softening, on the basis of the 
molecular field approximation, by taking into account 
the effect of spin-correlation and fluctuation. However, 
it is difficult to obtain the information of the excitation 
spectrum and the damping of the excitation under the 
high external field in the neutron scattering experiments. 
It is interesting to investigate the effect of softening of 
magnetic excitation and the spin dynamics in the order- 
ing process in the singlet ground state system. 

In these compounds, the magnetic nuclear spin-lattice 
relaxation time Ti is expected to be a useful probe for 
the study of spin dynamics of magnetic systems associ- 
ated with low lying excitations caused by the softening. 
Previously, the temperature and field dependences of Ti 
of ^^Rb have been measured in RFC.^^-* The effect of the 
soft mode was found in the temperature dependence of 
the relaxation rate T^^- Although the behavior of T^^ 
was explained qualitatively, quantitative understanding 
was not necessarily satisfactory. As for CFC, the field 
dependence of Ti of ^'^^Cs have been measured up to 
14 T {H//c) at the constant temperature of 4.2 K in the 
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paramagnetic region. The remarkable field dependence 

of was found with the maximum value at He in the 
field region where the magnetic phase transition occurs. 
The experimental data were interpreted from the view- 
point of the direct process, which occurs via the over- 
lapping of the ground state and excited state due to the 
intra-chain exchange interaction. However, the presence 
of relatively weak temperature dependence in Tj~^ has 
remained unsolved. 

The purpose of the present work is twofold: It is 
to clarify the relaxation mechanism associated with the 
softening, and then to obtain the information of the soft- 
ening of magnetic excitation in the paramagnetic state 
of the singlet ground state system. We have mcasiired 
the temperature dependence of Tj^^ of ^'^'^Cs in CFC, in 
the temperature range from 1.5 K to 22 K mainly ap- 
plied parallel to the c-axis up to 7 T. Preliminary results 
were already published. We have also measured the 
temperature dependence of Tj~^ of ®^Rb in RFC in the 
temperature range from 1.5 K to 22 K mainly applied 
perpendicular to the c-axis up to 3.3 T. The analysis 
is performed on the basis of the theory of DCEFA de- 
veloped by Suzuki. We shall survey the experimental 
procedures in §2 and results in §3. The derivation of the 
theoretical equation of T^^ is explained in §4, and §5 
is devoted on the discussions over the obtained results. 
Finally, the conclusion is summarized in §6. 

§2. Experimental procedures 

Single crystals of CFC (RFC) were made by mixing 
the equimolar materials of CsCl(RbCl) and FeCl2. The 
source materials were qualified for the purity of 99.9 %. 
Because of the hygroscopic nature of FeCl2, the materials 
were dehydrated using HCl gas in the ampoule. ""^^^ After 
the dehydration procedure, single crystals were grown 
by the vertical-gradient-freeze Bridgman method. The 
sealed ampoule was lowered by winding down a siltex 
code at a speed of 2cm/day. Single crystals as large as 
4^5 cm^ were obtained for the maximum size. 

The measurements of NMR were performed using the 
conventional phase-coherent pulsed-NMR spectrometer, 
which covers frequency range of 5 ~ 300 MHz. The 
NMR spectra were obtained by observing a spin echo 
signal as a function of an external field. The spin lattice 
relaxation time Ti was measured by observing a recovery 
of the signal intensity, after the saturation of nuclear 
magnetization by comb pulses. 

We determined the nuclear spin lattice relaxation time 
Ti of ^^Rb and ^^^Cs as follows. As for ^^Rb in RFC, 
the recovery curve of the nuclear magnetization (/ = 
3/2) is given as follows, when the quadrupole splitting is 
large; 

Mo - M(t) , t , , . 6t , ,^ , 

Mo = • + ' • •^"P^-tT^- 

Here the values of coefficients a and b were chosen as 
adjustable parameters, because it was difficult to ob- 
tain the complete saturation of magnetization. Reflect- 
ing the fact that the quadrupole splitting is sufficiently 
large compared to the spectral line width, the spectrum 
splits into three lines. Figure 1 shows a recovery curve of 
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Fig. 1. A recovery of the intensity of a spin echo signal of ^'^Rb 
in the paramagnetic state at T=2.8 K and H=3.33 T (HJ-c). 
The sohd hne is the fit to eq.(2.1) in the text. The inset shows 
the whole NMR spectrum of ^'''Rb. 



the intensity of a spin echo signal measured at the central 
line, and the whole NMR spectrum is shown in the inset. 
As is seen, the obtained curve is well fitted to eq.(2.1). 
On the other hand, in the case of CFC, the recovery 
curve is almost exponential as shown in Fig. 2. The inset 
in Fig. 2 represents the NMR spectrum of ^^^Cs. Since 
the quadrupole moment of ^"^^Cs is extremely small, one 
absorption line was observed as shown in the inset of 
Fig. 2, and the value of Ti is determined by fitting to the 
single-exponential recovery equation. 

§3. Experimental results 

Figure 3 shows the temperature dependences of the 
relaxation rate Tj"^ for CFC at various fields up to 7 T 
{H//c), and at 6.0 T (7J_Lc). As is seen, there is a dis- 
tinct difference in the temperature dependences of T-f^ 
between the lower (iJ<3 T) and higher {H>3 T) field 
regions. At low fields, T{~^ decreases remarkably with 
decreasing temperature. As the field increases, the tem- 
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Fig. 2. A recovery of the intensity of a spin echo signal of ^^^Cs 
in the paramagnetic state at T=20 K and H=7 T (H/fc). The 



inset shows the NMR spectrum of ^^^Cs. 
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Fig. 3. The temperature dependences of Tj ^ of ^^^Cs at various 
fields applied parallel (or perpendicular) to the c-axis. 



perature dependence of T-^^ becomes more moderate, 
and this suggests that the gap energy decreases. As is 
shown in the open circle in Fig. 3, the vahie of T^^ for 
HJ-C exhibits remarkable decrease with decreasing tem- 
perature reflecting the fact that the gap does not close. 
On the other hand, in the higher field region, T-^^ is 
almost temperature-independent down to 5 K, and ex- 
hibits an appreciable increase at lower temperatures. It 
seems that such a qualitative difference corresponds just 
to the absence and the realization of a magnetic phase 
transition. 

Figure 4 shows the temperature dependences of T^^ 
of ^^Rb in RFC. In this compound, a phase transition 
occurs under zero field at 2.6 K. As shown by the data 
obtained at 0.82 T {H//c), T^^ exhibits a drastic increase 
around the transition temperature. Above about 5 K, 
arc almost constant with the value 
In the case of H±c, the transition 



the values of 
of Tf 1 « 200 s- 
temperature decreases to zero K at 2.2 T with increas- 
ing external field. The relaxation rate measured at 0.75 
T and 1.53 T shows drastic increase around the phase 
transition temperature. It was found that T^^ increases 
gradually with increasing temperature above 5 K. While, 
as shown in the data taken at 3.33 T, T^^ for the field 
above the critical field shows drastic decrease with de- 
creasing temperature. 

As a common feature in both compounds, the relax- 
ation rate decreases remarkably with decreasing temper- 
ature in the field range where the phase transition docs 
not occur. On the other hand, in the field region where 
the phase transition occurs, the relaxation rate exhibits 



the divergent behavior. In the following, we shall ana- 
lyze the above experimental results on the basis of the 
theory of the DCEFA. 

§4. Theoretical background 

4-1 Description of electron spin system 

In the DCEFA theory, the spin product is decoupled 
in the following manner, 

Siij,Sjn — > Siij,{{Sj^) + a{Si^ — {Sin))) 

+ Sj^iiSin) + a{Sjn - {Sjn))), (m = X, y, z) (4.1) 

where (Sifj,) and (Sjfj,) are the thermal averages, and a 
represents a correlation parameter to be determined self- 
consistently. Due to the crystalline field, the electron 
spin system is largely influenced by the direction of H 
with respect to the anisotropy axis, the c-axis. In view 
of this, we treat individually the electron spin systems 
in the two cases of H//c and HLc. In the followings, we 
define the a;-axis and 2-axis as the crystalline a-axis and 
c-axis, respectively, in the hexagonal structure of CFC 
and RFC. 

4.2 H//C 

First we consider the case oi H//c. The effective single- 
ion Hamiltonian of the i-th site is expressed as"'^^-' 



H 



ESj. 



(4.2) 



with E ^ D + a(J^(0) - J/''(0)) and B^/ = g/fUBH + 
2J'''(0)(1 — a){Sz)- Here, {S^) is the thermal average of 
the induced moment which is independent of the z-th site; 
JT(0) = J7+3J2 (7 = -L, /) represents g = component 
of the Fourier-transform J'^{q) of the exchange integrals, 
where Ji and J2 are the intra- and inter-chain exchange 
interaction, respectively (we use the symbols /, _L for 
expressing the components parallel and perpendicular to 
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Fig. 4. The temperature dependences of Tj ^ of *''Rb at various 
fields applied parallel (or perpendicular) to the c-axis. 
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the c-axis). 

Figure 5 shows the energy level structure of the ef- 
fective single-ion Hamiltonian in eq.(4.2) for H //c. The 
ground state and the first excited state crosses at H^- 
The energy eigenstates of H^^ are given as |1) = |0), 



±> 



I0> 




Fig. 5. Schematic illustration of the energy level structure of the 
effective single-ion Hamiltonian for H //c given in cq.(4.2) in the 
text. Excitation energies rciHCScrited as Lo'q '' and ujq ' in the 
general case, originate primarily from the single-ion excitations 
|0) -> |±) for q = 0. 



|2) = |— ) and |3) = |+), and the corresponding eigen- 
states are given a.s Ei = Q, E2 = E — B and Ez = E + B, 
respectively (here, |0), |±) are the eigenstates of Siz). 
The matrix representations of (C = +> ~) with re- 
spect to the eigenstates, |0 >, |— > and |-|- > are given 
as follows; 



Si+ — 



Siz — 



V2 


V2 



) Si- — 



V2 
V2 





(4.3) 



It is noted that the magnetic excitation exists between 
the states connected by non zero value of the off-diagonal 
elements. According to the DCEFA theory, which takes 
into account the qt-dependence of magnetic excitation, 
the excitation energies defined as Wg"' and Wg"' in Fig. 5 
are obtained as poles of the dynamical susceptibilities 
X"^(q, w) = x'^'iQi'^)- The derivation of the dynamical 
susceptibilities x'^''(q,aj) {(,r] = -I-,— ,z) are shown in 
Appendix A. We note that x^''(q,cj) = x^^('Zi'^) origi- 
nates from the diagonal elements of the matrix represen- 
tation of Siz, and it is derived by considering the response 
of the induced moment to the infinitesimal displacement 
of the molecular field. 

In general, the relaxation rate Tj~^ is given by the 
Fourier-transform at the nuclear Larmor frequency uj^ 
of the fluctuating local field transverse to the nuclear 
quantization axis (defined as the Z-axis) as follows; ^''^ 



1 

7\ 



7^ 



dt{{h+{t)h- (0)}) exp(iw„f), (4.4) 



where {AB} denotes the symmetrized product \{AB + 
BA). The transverse local field at the "3Cs (S^Rb) nu- 



cleus is due to the dipole field coming from the surround- 
ing electron spins Sic^ of Fe^+ ions. Then the above equa- 
tion is rewritten in the following equation. 



1 



/OO 
, „ -CSO 



2N 

x({5^(i)5l„(0)})exp(ia;„i). 



(4.5) 



where A'^'^{q) is the geometrical factor which is calcu- 
lated as the product of the dipole tensor; Sq is the 
Fourier-transform of the electron spin operator Si(^ with 
respect to the wave vector q. 

As we described in the previous paragraph, we ob- 
tain the eigenfrequency Wg' (^ = -|-, — ) of the singlet 
ground state system. Here, we extract explicitly from 
the time dependence of the spin operator the part re- 
lated to We use the fluctuation-dissipation theorem 
between the spin correlation function and the relaxation 
function as;^*-* 

/OO 
dte'^\{5S<q{t)5Slqm) = Ep{u;)R'i^{uj), (4.6) 
-00 



with 



£;,(u.) = -coth(^) 

,(«) 



(4.7) 



where w is taken to be ujq or uj^ according to the relevant 
components, and Rqito) is the Fourier transform of the 
relaxation function Rqit). We here assume that Rq(t) 
in eq.(4.6) is described in terms of q-dependent static 
susceptibility x'"'' (q) and normalized relaxation function 
fq\t) as 



(4.8) 



The static q-dependent susceptibilities x''''(q) are evalu- 
ated by putting w = in the dynamical susceptibilities 
X''''(<7, tij). Then, the dynamical nature in x''''('7:'^) is 
transferred to the normalized relaxation function, and 
the problem is reduced to the evaluation of fq [t) . 

Now we look at the normalized relaxation function in 
eq.(4.8). As shown in the previous paragraph, the ex- 
cited state of singlet ground state is charactrized by the 
excitation energy ijjq\ Then, when the damping effect is 
ncglectcid, the time evolution of the excited state is de- 
scribed as cxp(ici;q't). Let us take into account the damp- 
ing effect which is expressed by the exponential function. 

As for the fluctuation transverse to the external field, 
wc define fq{i) and disregard the effect of the Tipper 



branch 



J- 



Treating the time evolution of the lower 



branch ujq' which is effective to the nuclear magnetic 
relaxation, the normalized relaxation function may be 
given as follows; 



/^(t) = exp(^c^<-'^-^^^). 



(4.9) 



Physical meaning of Fq can be considered to be the 
lifetime of magnetic excitation ijjq \ and broadening of 
the spectrum occurs due to Fq. The concept of life- 
time is not considered within the frame of the DCEFA 
theory; however, when there is interaction between the 
excitations, or the higher order terms of interactions are 
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included, then the magnetic excitation may have Hfetime 
whose value is given by the inverse of Fq . 

As for fq{t), considering that this is the fluctuation 
of the induced moment along the external field, we may 
put 



exp{-r^t). 



(4.10) 



Arrangement of the induced moment might be disordered 
by the propagation of magnetic excitation (or the trans- 
verse fluctuation), then this could be the physical origin 

ofrl 

Here we define /q(w) as the Fourier-transform of fq{t) 
with respect to w. Figure 6 shows an illustration of the 
frequency dependence of fq{uj) and ,fq(uj) near uj = 0. 
fq{uj) and /q (w) have their centers at zero and Wg"', and 
broadening of the spectra are characterized by Fq and 
Fq. As shown in Fig. 6, the relaxation rate will be ex- 




Fig. 6. An illustration of the frequency dependence of fq (u)) and 



plained in terms of two different kinds of electron spin 
fluctuations, and the ratio of each contribution will be 
determined by the values of Fq and Fq (the effect of the 
excitation of Wg"' mode can be ignored). If there is no 
broadening of spectrum due to Fq, then fq{uj) do not 
have frequency of Un component, and this fluctuation 
will not work as a thermal reservoir of the nuclear spin 
relaxation. 

Referring to the above considerations, the theoretical 
equation of Tf^ is given as Tf ^ = (T^^)^^ + (Tf ^)-^ with 

q 



{Try = 



2 N 



(4.11) 



(7^r')" = |^E(5-MB)M-(9) 



2 N 



X coth( 



fill) 
~2 



flLO 

2KT 



)x"(q) 



2r^ 



(4.12) 



where Fq 
UJq'' + U)„ 



^ ojn is assumed for eq.(4.11), and lo = 
in eq.(4.12); considering that Uwn <^ k^T 
for the relevant temperatures, we put Ep{ijj-a) ~ k^^T 
for eq.(4.11). The expressions of the geometrical factors 
A~'{q) are shown in Appendix B. 



i.3 H±c 

Next we consider the case for H //x (crystalline a-axis), 
which corresponds to the present experimental condition. 
The effective single-ion Hamiltonian is expressed as^^^ 



(4.13) 



with = q^hbH + 2J-"(0)(1 - a){S^). Here {S^) is 
the thermal average of the induced moment. 

Figure 7 shows the energy level structure of the effec- 
tive single-ion Hamiltonian. The gap energy between the 



±> 
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Fig. 7. Schematic illustration of the energy level structure of the 

effective single-ion Hamiltonian for H^c {H //x). Excitation en- 
ergies represented as a;^ ' , u)q ' and ' in the general case, origi- 
nate primarily from the single-ion excitations |1) |2), |2) — > |3) 
and |1) |3) for q = 0. 



ground state and the flrst excited state increases with 
increasing external field. The eigenstates of H^^ , which 
are given by the linear combinations of |±) and |0), and 
the corresponding energy eigenvalues are given as fol- 
lows, 

|1) = di/V2(|+) + h)) + V2ci|0); = 1/2{E - W), 

(4.14) 

|2) = (1+) - |-))/%/2; E2 = E, (4.15) 

|3)=c,(|+) + |-))-4|0); E3 = l/2{E+W), (4.16) 

where W = {E'^ + AB'^f/^, F = {{E + Wf + 4B^y/^, 
c, = V2B/F, d, = {E + W)/F, C2 = djV2 and 4 = 
•\/2ci. The matrix representations of Si^ {fi = x,y,z) 
with respect to the eigenstates, |1), |2) and |3) are given 
as follows; 



2A/2diCi V2{d^c^ - d^c, 



-2\/2diCi 



ly 











V2(di 


Ca — d2 


ci) 





idi 


' \ 


—idi 





id2 





-id2 


/ 





V2c, 










V2c2 





\/2c. 






(4.17) 



According to the DCEFA theory, excitation energies are 



Uq , ujq 

as poles of x^*'(g,w). 



,(3). , ,(1) 

jq , ujq 



and 



are obtained 



X^''(9,w) and x^^(g,t^), and 
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is obtained as the pole of x^^(q,w). The expressions of 

the dynamical susceptibilities x^'^iQ:^) (p-t'^ = 2;,?/, 
are given in Appendix A. The terms of x^^ili^) 
X^^(q,w) can be expressed as a sum of two terms which 
have poles at tjg' and a;^' respectively, as shown in 
eq.(A-13). As for x^^il^ i^)) two expressions are given for 
(jf^{u>) according to the case of w = and 7^ as shown 
in eq.(A-6). The former case results from the diagonal 
elements of Six in eq.(4.17); we define xZaiQ^ which is 
derived by substituting the expression of {ui = 0) into 
eq.(A-7). The latter case results from the off-diagonal 
elements of Six in eq.(4.17), and tjg' is obtained as the 
pole of which is derived by substituting the 

expression of 4)^^^{oj ^ 0) into eq.(A-7). 

As for the normalized relaxation function, we use the 
relation in eq.(A-13) and treat the time evolution of the 
transverse spin fluctuation separately with respect to a;^ ' 
and ojq\ Then the normalized relaxation function for 
the dynamical susceptibilities Xw"{l^^) Xu^^Qi^) 
{n, v = y,z) may be given as follows; 

f^'\t)=eMi^'^'t-r^'% (4-18) 
f^\t)=eMi'^qt-r'qh). (4.19) 

As for the fluctuation x^^{Qi ^) which is along the ex- 
ternal field, we do not take into account of the effect of 
the ' mode; because it does not show softening and the 
excitation energy increases largely by applying an exter- 
nal field. On the other hand, we take into account of 
the longitudinal fluctuation of the induced moment, and 
the normalized relaxation function for this fluctuation is 
defined as 

/^-(f)=exp(-r^'t), (4.20) 

where suffix (0) is used to distinguish the damping con- 
stant of the longitudinal fluctuation for H±c from that 

of H//C. 

Then the theoretical equation of Tj~^ is given as T^^ = 
(Tf i)i + (Tf 1), + (Tf i)o with 

(^r')^ = YJ^E^<^otH^){{g^i,^rAyy{q)xl\{q) 
q 

2j-(i) 

+ ig,^,fA-{q)x:{{q)}j^;^j^^^^, (4.21) 

q 

2/^(2) 

+ (S//MB)'^^^(g)xr.(9)}--I5u^^— ^, (4.22) 

(rq'r + {coy 

(Ti')o = ^^E^-^(»-'^-)'-^''(9)^-S(9)7i> (4-23) 

where w = w^''-fwn for eq.(4.21), w = t^g'-l-tJn for eq.(4.22) 
and Tq"' ^ uj^ is assumed for eq.(4.23). The expres- 
sions of A^'^{q) are shown in Appendix B. From the re- 
lation Ay^{q) = A'^yiq), and x^'iQ,'^) = -X^^(9,w) in 



eq.(A-9), we obtain the following equation, 

Ay^iq) ■ x^'{q,^)+A'y{q) ■ X^^(g,a;) = 0, (4.24) 

thus X^^(Q)'^) and X^^(9)^) do not have any contribu- 
tion to the relaxation rate. 

§5. Analysis and Discussion 

In the previous section, the expressions of the relax- 
ation rate have been shown as the sum of the contri- 
bution of diff'erent magnetic modes. We introduced the 
damping constants for each magnetic excitation in the 
singlet ground state system. It is diflicult to evaluate 
these factors, theoretically; then, we determine the val- 
ues of damping constants by the fitting of the theoretical 
curve to the experimental results of T^^. The numeri- 
cal values of Hamiltonian parameters for CFC and RFC 
are shown in Tabel I, which have already been evaluated 
by Suzuki and Makino from the experimental results of 
ESR and neutron-scattering. 



Table I. The Hamiltonian parameters of CFC and RFC evalu- 
ated by Suzuki and MaJsino.^^^ 





CsFeCls 


RbFeCls 


9// 


2.54 


2.54 


9± 


3.84 


3.84 


D/k^ (K) 


20.1 


21.3 


J^/k^ (K) 


3.72 


5.05 


J.j-/kB (K) 


-0.201 


-0.429 


\-h/-h\ 


0.055 


0.085 


J /./- 
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5.1 H//C 

First, we consider the case of H//c in CFC. We show 




Fig. 8. (a) Temperature dependence of x^^ io) along q = 
(1/3, 1/3, Qc) and (b) X''~{q) along q = {q,q,0), calculated for 
CFC at if=5.0 T on the basis of the DCEFA. 



typical example of the q dependence of the static sus- 
ceptibility x''{q) in Fig-8, calculated at iJ=5.0 T. Figure 
8(a) represents typical example of the q dependence of 
X^^{q) along q = {1/3,1/3, Qc), each components being 
along the reciprocal-lattice vectors a*, b* and c*, re- 
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spectively. As is seen, x''{q) decreases appreciably with 

decreasing temperature, which resuhs from the behavior 
of the longitudinal single-ion susceptibility (p" {lo = 0) 
characterizing the singlet ground state system. Figure 
8(b) represents typical example of the q dependence of 
X^(q) along q = {q,q,0). As is shown, X'^{q) exhibits 
a distinguished tendency of divergence around the K- 
point with decreasing temperature, which is associated 
with the softening of magnetic excitation oj'q ' . 

As an example, we show the q-depcndcncc of A^' {q) 
along q = (1/3, 1/3, grc) and q = {1/2, 1/2, Qc) in 

(a) (b) ^ 




0.1 0.2 0.3 0.4 0.5 0.1 0.2 0.3 0.4 0.5 



Fig. 9. (a) Geometrical factor y.AII (q) along q = (1/3, 1/3, Qc) 
and q = (1/2, 1/2, Qc), and (b) r'^ x A^{q) along q = {q,q,0) 
where r is the distance between the ^^^Cs nucleus and the near- 
est Fe^+ ions (r = 4.4 x 10~® cm). 




5 10 15 

r(K) 



Fig. 10. The results of fitting at H=1.Q T and 2.2 T for CFC. 
The solid lines and the dotted lines represent the fitted curves 
obtained from eqs.(4.11) and (4.12) in the text. The heavy lines 
represent the total relaxation rate = (Tf^)// + (T~^)^. 



Fig.9(a), and A^{q) along q = {q,q,0) in Fig. 9(b). 

It should be noted that A'^ (q) is zero and A^(q) takes 
the maximum at the K-point, which corresponds to the 
120° spin-structure in the c-plane of the ordered state. 

In the following, we analyze the experimental results 
of Tf 1 in CFC in terms of the contribution of (T^^)^^ 
and (Tj"^)-^ (eqs.(4.11) and (4.12)), with two variable 
parameters Fq and Fq . Here we assume that F^ IS m- 
dependent of q and temperature, then the symbol Fq is 
expressed as f^'^ . As for Fq, we assume that Fq is in- 
dependent of temperature, and considering that x^(iir) 
increases largely and exhibits divergent behavior at low 
temperatures at the if -point, we put Fq fa F^. 

For the comparison, we show the analysis for the re- 
sults in the low field region and high field region, re- 
spectively. First we look at the low field data which 
are shown in Fig. 10. The heavy lines in Fig. 10 repre- 
sent the best fitted curve of the total relaxation rate 
Tj^^ = (Tf ^)-'/ + (Tf ^)-^. The solid lines represent the 
fitted curves of (T^^)^^ obtained by choosing a constant 
value of UF'' /k^ = 1.3 K, so as to fit the experimental re- 
sults at low fields. Hereafter we leave out the factor Ti/k^. 
As is seen, the agreement between the data and the fitted 
curves of (Tj~ '^)" is satisfactory over the whole tempera- 
ture range at 1.0 T. Next, we look at the high field results 
which are shown in Fig.ll. The solid lines represent the 
calculated curves of {F^^)'' which are obtained by using 
the same value of F" — 1.3 K. As is seen, the discrepancy 
between the experimental results and the calculation is 
serious. It is difficult to explain such a discrepancy by 
adjusting the value of F'' , because the diverging behav- 



CsinCsFeCl - 
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Fig. 11. The results of fitting at if=3.5, 4.0, 5.0 and 7.0 T for 

CFC. The solid lines and the dotted lines represent the fitted 
curves obtained from eqs. (4. 11) and (4.12) in the text. The heavy 
lines represent the total relaxation rate T^^ = (T^^)// -|-(Tj~^)-'-. 
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ior in Tf at low temperatures cannot be obtained from 

the calculatioon of (T^^)^^ . Such a discrepancy may be 
associated with another term of {T^^)^. The dotted 
lines in Figs. 10 and 11 represent the theoretical curves 
for (T-|~^)^, which are obtained by choosing the values 
of so as to fit the experimental results at low tem- 
peratures. The heavy lines in Fig. 11 represent the total 
relaxation rate Tf ^ = (Ti^)" + (^f ^)^. Then we find 
that the reasonable fit is attained by taking that in- 
creases with increasing an external field, while a constant 
value is expected for F^' . The best fit values of Fj^ at 
each field are summarized as follows; Fj^ w 0.02, 0.08, 
0.4, 0.6, 1.6 and 2.0 K (for 1.0, 2.2, 3.5, 4.0, 5.0 and 
7.0 T). It turns out that there is a significant difference 
in the value of Fj^, which is very small in the low field 
region, and relatively large in the high field region. This 
fact suggests that the damping of magnetic excitation 
becomes remarkable in the high field region where the 
effect of the softening is more pronounced as the sign of 
the occurrence of 3D-LR0. 

Next we consider the results of Tj~^ in RFC. The anal- 
ysis was pcirformed following the same procedure in CFC. 
The value of F^'' is determined so that the contribution 
of {Tj^)^^ explains the data at higher temperatures. The 
contribution of (F^^)^ is determined so as to fit the 
curve of the total relaxation rate to the experimental 
results around Tn- At ff=0.82 T, the best fitted curves 
are obtained for F^/ =3.0 K and F^ =0.4 K. As shown 
in Fig. 12, the agreement between the theoretical curves 
and the experimental results is very good. 

5.2 H±c 

In the case of HJlc, we first consider the results in 
RFC. It is noted that the effect of the external field is 
more evident than in the case of CFC. We show the typ- 




Fig. 12. The fitted curves ofTj"^ at H=0.82T for RFC. The solid 
Une and the dotted line represent the fitted curves obtained from 
eqs.(4.11) and (4.12) in the text. The heavy line represents the 
total relaxation rate Tf^ = {T~^)" + {T~^y . 



(a) (b) 




0.1 0.2 0.3 0.4 0.5 1 0-2 0.3 0.4 0.5 

1, 9= (q^.q^.o) ^=(q,.qe'0> 



Fig. 13. (a) Temperature dependences of x!^{q) along q = 
{q, q, 0) and (b) X^^iq) along q = (q, q, 0), calculated for RFC at 
H=1.0 T on the basis of the DCEFA. 



ical example of the q dependence of xZoil) ^^'^ X^^{q) 
along q = {q,q,0) in Fig. 13, calculated at H=1.0 T. 
As is seen, xZ^{q) dcc;rcases appreciably with decreas- 
ing temperature, which mainly results from the behav- 
ior of (P^^{lu — 0). The temperature dependence of 
(Tj~^)(, follows the behavior of xZoil)'^^ ^^'^ decreases 
with decreasing temperature. On the other hand, as seen 
in Fig. 13(b), x^^{q) increases with decreasing tempera- 
ture, and exhibits a distinguished tendency of divergence 
around the K-po'mt; this behavior is associated with the 
softening of magnetic excitation ujj^ and mainly comes 
from the contribution of x^^(K'). 

Figure 14 shows the calculated results of ^^^(q), 
Ayyiq) and yl^^(q) along q = {q,q,0) (Fig.l4(a)) and 
q = (1/3, 1/3, gc) (Fig.l4(b)). It is noted that A^^iK) 
has large value around the i^-point; therefore, the ef- 
fect of the divergence of xUi{K) will largely appear on 




0.1 0.2 0.3 0.4 0.5 0.1 0.2 0.3 0.4 0.5 



^= 9 = (q^,q^,0) <!r = (1/3, l/3,qp 



Fig. 14. (a) Geometrical factors x A'^''{q) along q = {q,q,0) 
and (b) X ^'"'(q) along q = (1/3, 1/3, gc), where r is the 
distance between the *^Rb nucleus and the nearest Fe'^^ ions 
(r = 4.3 X 10-^cm). 
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In the following, wc analyze the experimental results of 
Tf 1 in terms of (T^^),, {T{% and {T^\ (cqs.(4.21), 
(4.22) and (4.23)), with T^"', T^^' and r^"> as an ad- 
justable parameters. We here assume that r^°' and T^"' 
are independent of q and temperature, thus the symbols 
r^°' and r^'> are expressed as r<°' and r<=^). As for T^", 
we assume that Pg ' is independent of temperature, and 
considering that xZ\{I^) increases largely and exhibits 
divergent behavior at low temperatures at the ii'-point, 

we put r^'' « r^j^. 

First, wc assume that the contribution of {T^^)2 is 
very small, because this is due to the excitation between 
the excited states. Figure 15 shows the fitted curves for 
the experimental results at 1.53 T and 3.33 T.When the 
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Fig. 15. The results of fitting at //=1.53 T and 3.33 T {H±c) 
for RFC. The sohd-Unes, the dashed-dotted lines and the dotted 
lines represent the fitted curves obtained from eqs.(4.21), (4.22) 
and (4.23) in the text, respectively. The heavy lines represent 
the total relaxation rate T^^ = (T^^)i + (Tf + (T-^'-)o- 



when the value of Fj^ is varied from 0.4 K to 1.5 K, and 

we take it that damping of mode becomes severe as 
the sign of the occurrence of 3D-LR0. 

Now we turn to the data at 0.75 T, and fit the theo- 
retical curves using the values of F'"' and 7^^' obtained 
at 1.53 T. However, it turns out that the reasonable fit 
cannot be obtained at 0.75 T within the contributions 
of (Tj7^)i and {T^^)a, especially at higher temperatures. 
Considering that Wq' is the excitation between the ex- 
cited states, and {T^^)2 increases with increasing tem- 
perature, now we take into account of the contribution 
of {T,-%. 

Figure 16 shows the fitted curves for the experimental 
results at 0.75 T. The dashed-dotted line represent the 
fitted curve of {T^^)2, and the best fit is obtained when 
the contribution of (Tf^)^ is included with T*"' = 0.04 K. 
It is conjectured that magnetic excitation of Uq^ mode, 
which do not exist in case of H//c, would have appreciable 
effect on T^^. In Fig. 15, we have also shown the fitted 
curves of {T^^)2 and the total relaxation rate at 1.53 T 
and 3.33 T with F^^^ = 0.04 K. 

Next wc analyze the experimental results of Tj~^ of 
^^^Cs in CFC. Figure 17 shows the fitted curves of T^^ 
at 6.0 T. Referring to the analysis for the data at 3.33 
T in RFC, the analysis was performed comparing the 
contributions of (Tj~^)o and (Tj~^)i, because the contri- 
bution of (T'j~^)2 can be still more disregarded owing to 
the larger effect of the external field. The value of T*"' is 
determined so that the contribution of (rj~^)o explains 
the data at higher temperatures, and the value of F'^^ is 
determined so as to fit the curve of the total relaxation 
rate to the data around 2 K. The value of T*^' follows 
that obtained in RFC, because this value will not be 
largely different between these materials. The best fit- 
ted curves are obtained with F^^ = 0.07 K, T'"' = 0.04 
K and T'"' = 2.5 K, respectively. 



external field of 3.33 T is applied in the c-plane, the ex- 
citation energy increases largely, and we may ignore the 
contribution of {T^^)^. The solid line and the dotted 
line represent the fitted curves of (Tj"^)i and {F^^)o, 
respectively. The value of {F^^)o decreases, while that 
of (Tj~^)i increases with decreasing temperature; there- 
fore, at 3.33 T, the analysis was performed assuming 
that the relaxation rate is mainly due to the contribu- 
tion of (rj~^)o at higher temperatures, and the descrep- 
ancy from the data at low temperatures is ascribed to 
the contribution of {T^^)^. The best fit for the data at 
3.33 T is obtained with r<°' = 5.0 K and T]^' = 0.4 
K. The heavy line represents the total relaxation rate 
Tj-i = (Tf ^)i + (Tf ^)2 + (rf^)o, and the estimation 
of the small contribution of (Tj~^)2 will be explained 
later. At 1.53 T, the phase transition occurs, and the 
tendency of divergence of T{~^ appears at around Tn- 
We expect that this behavior results from the contribu- 
tion of {F^^)i. Refering to the fitting procedure of F^^ 
in CFC, we expect that the damping constant of the soft 
mode {F^) depends on the field and the occurrence of 
the magnetic phase transition. The best fit is obtained 




riK) 



Fig. 16. The results of fitting at H=0.75 T {H±c) for RFC. The 
solid-line, the dashed-dotted line and the dotted line represent 
the fitted curves obtained from eqs.(4.21), (4.22) and (4.23) in 
the text, respectively. The heavy line represents the total relax- 



ation rate Tf ^ = (Tf i)i + (Tf 1)2 + (Tf ^)o. 
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5.3 Discussion 

Throughout the analyses based on the DCEFA, we 
have evaluated quantitatively the damping constants of 
the spin fluctuations which have effects on the relax- 
ation mechanism. The damping constants for the fluc- 
tuation of the induced moment are determined at some- 
what smaller values than those of the relative exchange 
interactions. The lifetime of magnetic excitation is long 
and the values of damping constants are evaluated more 
or less about 0.1 K, in the fleld region where the sys- 
tem remains non-magnetic down to zero K. On the other 
hand, in the fleld region where the ordering occurs, it was 
deduced that the lifetime of soft mode is short, and the 
values of damping constants are estimated as large as 1.5 

2.0 K. If we consider that there is interaction between 
magnetic excitations, then the collisions between exci- 
tations will explain the field dependence of Fj^ (^jc)- 
When the magnetic phase transition occurs at low tem- 
peratures, then the times of collisions might increase at 
T ~ Tn; because large numbers of the excitations over- 
come the activation energy when the excitation energy 
decreases to zero. Then owing to the effect of the an- 
harmonic terms of the interactions, lifetime of magnetic 
excitation will become short and the value of damping 
constant will increase. 

The evaluated values of damping constants are sum- 
marized in Table II. It is very instructive to refer to the 
experimental results of inelastic neutron scattering on 
CFC, which have measured the line width of the excita- 
tion spectra under various external fields applied parallel 
to the c-axis at r=1.6 K.^°) The line widths of mode 
increases largely in the field region where the magnetic 
phase transition occurs, and the absolute value of the 
line width is about 1 K above 4.5 T. This value is only 
different from the estimated value of for the factor 
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Fig. 17. The results of fitting at _H"=6.0 T (HLc) for CFC. The 
solid-Hne, the dashed-dotted fine and the dotted fine represent 
the fitted curves obtained from eqs.(4.21), (4.22) and (4.23) in 
the text, respectively. The heavy line represents the total relax- 



of 2 at the maximum. 

We emphasize that various values of damping con- 
stants shown in Table II are new informations, which 
are obtained for the first time through the analyses of 
Tfi in the singlet ground state system. In particular, it 
is noteworthy that we could estimate the damping of the 
longitudinal fluctuation and the excitation oiw'q modes, 
and show their effects on T^^; inelastic neutron scatter- 
ing technique is difficult for investigating the fluctuation 
at w Ri 0, and the excitation between the excited states. 
Informations of such fluctuations are obtained through 
the advantages of NMR experiment. 

§6. Conclusion 

In conclusion, we have measured the temperature de- 
pendences of T]~^ of ^^^Cs (^^Rb) in the singlet ground 
state system CFC (RFC) under various fields (i?<7 T) 
applied parallel or perpendicular to the c-axis. In the 
field region where the system remains non-magnetic 
down to zero K, it was found that a remarkable de- 
crease of Tj"^ was governed by the longitudinal fluctu- 
ation along the external fleld, which depends strongly 
on the appreciable energy gap between the ground state 
and the excited states. On the contrary, in the fleld 
region where a magnetic ordering occurs, an apprecia- 
ble increase of below about 5 K is ascribed to the 
transverse fluctuation perpendicular to the external field, 
which is associated with the softening of magnetic exci- 
tation. The damping constants of the normalized relax- 
ation function were evaluated for the fluctuations which 
have signiflcant effects on the nuclear magnetic relax- 
ation. It was deduced that the damping constant of the 
soft mode increases largely in the field region where the 
softening occurs, and then the tendency of divergence of 
T^^ is ascribed to that of the static susceptibility. In 
RFC, it was conjectured that the magnetic excitation of 
uj'q mode would have significant effects on Tj~^ in case 
of HLc. 

We thank Mr. S. Ueda for cooperation with the prepa- 
ration of sample. 

Appendix A: Dynamical Susceptibitlies 

The dynamical single-ion susceptibilities are calcu- 
lated by the following expressions; 



{Pn - Pm) < m\Sifj,\n >< n\Si,y\m > 



^ + Em — En 



(A-l) 



Table II. 
fit. 



The values of damping constants evaluated through the 





CsFcCls 


RbPeCls 


rk (K) 


0.02~2.0 


0.4 


r« (K) 


1.3 


3.0 


^ic'' (K) 


0.07 


0.4~1.5 


r(2) (K) 


0.04 


0.04 


r(o) (K) 


2.5 


5.0 



ation rate Tf ^ = (Tf i)i + {T-'^)2 + (Tf i)o. 



with pn = exp{-En/kBT)/Y^^exp{-Em/kBT) {n,m 
1 . 2 and 3). For the diagonal elements n = m, which 
exist for Siz in case of H//c and Six in case of ff_Lc, 
the dynamical single-ion susceptibilities are derived by 
considering the effect of the infinitesimal displacement of 
the effective field. Then for the case of H//c, wc obtain 
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k^T 



)Su.,o, (A-2) 



and the dynamical susceptbilities w) and x-^(q,a;) 

for the whole system are expressed in terms of the single- 
ion susceptibilities as follows: 



1 



(A-3) 



■p7(j7(q)-aj7(0))<^7(a;)' 

where p'' = 2 or = 1. For example, the expression for 
X^(q,w) is rewritten in the following form; 



(a; + c.^-')(o.-o.<y+') 



(A-4) 



Pn 1 -^mn 



En and w, 



(±) 



IS 



where Pmn = Pr, 

expressed as 

a;^±'= [E^ - 2E{J^{q) - aJ-(0))(2 - 3(5,)^) 
+{Sz)\J^{q) - aJ-(0))2](V2) ± lg>/^^H 
+2J//(0)(1 - - (J-(q) - aJ^(0))(5,)] . 

(A-5) 

In the case of HJ-C, the matrix representation of 

Six have both off-diagonal and diagonal elements, then 
the dynamical single-ion susceptibilities are respectively 
given as 

2p,,E,,iE/W)^ 

io-'-Ei 



+ 



2p,,{E/Wf 

^ [{2BIWf{p,+p,)-{Sx)\{i^ = 0). 



4^^(00), <j>yy{uj) and are derived by using eq.(4.17) 

and cq.(A-l), and other components of the single-ion 
susceptibility are vanishing. The dynamical susceptbili- 
ties x'^'^iQji^) for the whole spin system are expressed in 
terms of the single-ion susceptibilities as follows: 



+ 2J^{q)r'{co)cl>^y{iv)]/D{q,iv), (A-S) 
X^'{q, = cl>y^{ij)/D{q, to) = -x'y{q, u), (A-9) 

X'^iq,^) = - 2J^{q)cf>yy{u,)) 

+ 2J-(9)0^^(a;)</.^^(a;)]/D(q,w),(A-lO) 



(A-11) 



JZ{q) = JHq)-ar{0), 
Diq,uj) = 1 - 2J^{q)<j>yyiu;) - 2J^{q)r^{u;) 

+ AJi{q)J^{q){4>yy{u;W^{uj) - </.^^(a;)<^^^(a;)). 

(A-12) 

For example, x'"'(9,w) can be expressed as follows, 



x'Ji{q,oj) + xlliq,oj) 

jyy 



+ 



tVV 

-'2 



where ^ and are given as 



(A-13) 



n 



vv 



K')^-«)^' ' 



and ai, 61 are given as 

fli = 2{p,,(f,E,, + p^,dlE^^) - %pl,cldl2J" [q) 



(A-M) 



(A-15) 



K = ElEl2{'^ + 

-£^21 



E^. 



) 



X [1 -4(^ + ^)2 J// (q)]. (A.16) 



-£^21 



E,. 



,(1) 



and 



,(2) 



(jjq are given as 



,(1) 



(A-17) 



where Ai and Bi are given as 

Ai = + El) - {p.^E.^dl + p,,E,,dl)2J\q) 



(A.6) - {P..E,,(?, + p,,E,,(?,)AJ''{q) + lQpl,(?,dlJ'' {q)J^{q 



(A-18) 



B,=ElEl[l-2{ 



Pl2df _^ P23dl 

E,, 



)2J^{q) 



X [1- 4(^ + ^)2 J//(g)]. (A.19) 

Appendix B: Geometrical Factors 

The calculations for the geometrical factors were made 
by taking the dipolar sum only over the nearest-neighbor 
sites of the Fe^+ ions, considering the fact that this factor 
decreases as inverse sixth power of the distance between 
the ^■''■''Cs (^■^Rb) nucleus and the Fe2+ ion. 

For the convenient, wc consider the geometrical figure 
as shown in Fig. 18. The position of ^■'■^Cs (^''Rb) is taken 
at the origin, and the sites of the Fc^+ ions are placed 
on the vertexes of the trigonal prism, which are labeled 
from A to F. The coordinates of each site are given using 
the lattice parameters a and c. The /i'-component of the 
dipole field (hi) at the ^^^Cs (*''Rb) nucleus, due to the 
/i-component of the magnetic moment (mj) of Fe^"*" ion 
at the site i, is given as follows; 

K =i-r- S-^X = rf^'^rOmf , (B-1) 



P23C2 
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Fig. 18. The illustration of geometrical figure of the -'^'''^Cs (*''Rb) 
nucleus and the nearest neighbor Fe'^+ ions. The position of the 
^^^Cs (^'^Rb) nucleus is taken at the origin, then the sites of the 
Fe^"*" ions on the vertexes of the trigonal prism are labeled from 
A to F. The coordinates of each site are given using the crystal 
parameters a and c. 



where Vi is the distance from the origin to m,, and 

dfj,'p,{ri) is the dipolo tensor. Then the correlation of the 
fluctuating local field is related to the spin correlation 
via the product of d^'^(ri) as 

(J/i^'(i)(5/i'^'(0)) = J2 d^,^{u)d,.,{r,){5m'i{t)5m)m. 

(B-2) 

The geometrical factor A^'^{q) can be calculated as the 
Fourier transform of d^i fj_{ri)d^i,^{rj) with respect to the 
wave vector q. The nuclear Zeeman axis {Z) coincides 
with the .2-axis for H//c, therefore, Sh^ in eq.(4.4) is 
defined as Sh"^ = Sh^ ± iSh^. 

In the case of H//c, the expression of A^{q) is given 

as 

-4" (9) = {(i - + ^} ■ {6 + 6cos(27rg,)} 



4r, 



r}-{4cos(27rga + 27rg'6) 



A'^^{q) and A^^{q) are given as 
+ (t^To - ^TTo ) ■ {4cos(27rgc) - 2cos(27rg'a + 27rgfe 



+ 2 cos(27rga + 27r<75 — 27r(jc) + 2 cos(27r<7a + 2TTqb + inqc) 
+ 4cos(27rga) + 2cos(27r(ja — 2TTqc) + 2cos{2TTqa + 2-Kqc) 
+ 4cos(27rg6) + 2cos{2iTqb — 2iTqc) + 2cos(27rga + 27r(7c)}, 

(B-3) 

where q = {qa^lbilc)] c are the lattice parameters 
shown in Fig.18 and = (a^/3 + c'^/iy/'^. A/^{q) is 
given as 



A^q) 



16r, 



10 



{6—6 cos(27rg'c) — 2 cos(27rga + 27rg'b) 



+ cos(27r<7a + 2T:qb — 2'Kqc) + cos(27r<7a + 2'Kqb + 2'Kqc) 

— 2 cos(27rga) + cos(27rga — 2'Kqc) + cos(27rga + 2'Kqc) 

— 2cos(27rg(,) + cos{2nqb — 2nqc) + cos(27rQ'a + 27rgc)}- 

(B-4) 

In the case of 7?_Lc, the nuclear Zeeman axis (Z) coin- 
cides with the a;-axis, therefore, 5h^ in eq.(4.4) is defined 
as Sh^ = Sh^ ± iSh^. Then the expressions of A'^'^{q), 



2'Kqc) - 2 cos(27rga + 27rqfa - 2'Kqc)}, 



(B-5) 



1 2 o 2 2 

A''{Q) = {{^^ - - ^} ■ {4cos(2.,,) 



+ 2 cos(27rga + 2'Kqb - 2'Kqc) + 2 cos(27rg'a + 27rg6 + 27rgc)} 
+ {(i - + ^} • {4 + 4 cos(27rg„ + 27rg6)} 

,,1 a^.o a^c^ . ^ ,,1 a^.o a^c^ . 

.2cos(2.,.) + {(-3-^).(-3--,)-^} 

• {4cos(2^Q,) + 4cos(2^<jfc)} + {(1 - " ) . - " ) 

'0 



32r, 



10 



} • {2 cos(27rga - 2'Kqc) + 2 cos(27rga + 2'Kqc 



+ 2 cos(27rg6 — 27rgc) + 2 cos(2'7rg'6 + 2'Kqc)}, 



(B-6) 



+ 2 cos(27r(jf„ + 27rgb - 27rgc) + 2 008(277^0 + 27rg6 + 27rgc)} 



' 



2cos(27rgc) + {(^-^) 



1 _ Sc^ 2 _ 3a^c2 



} • {4cos(27rga) 



+ 4cos(2.5.)} + {(^-^)^ + |^} 

• {2 cos(27r(7a — 2Trqc) + 2 cos{2Kqa + 2'Kqc) 
+ 2 cos{2nqb - 2'Kqc) + 2 cos(27rqb + 2Kqc)}. (B-7) 
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